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I. INTRODUCTION 



This work is an extension of some analytical results that have already been presented in 
a previous paper on Brownian motion (Gaioli et al, 1997), hereafter referred as paper 
In I we have considered a model consisting of a Brownian oscillator of frequency Q linearly 



Equations numbered as {jf) in paper I will be labeled as (!•#)• 
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coupled to a bath of harmonic oscillators in the rotating wave approximation. This model 
has an exact solution from which we have studied the time evolution of the relevant physical 
quantities, e.g. the mean position (X(t)) and mean population (Na(t)) of the Brownian 
oscillator. We have shown that the equation of motion governing the evolution of (X) 
is a generalized- local-in-time form of the Langevin equation (in mean values), with time- 
dependent coefficients (see also Garcia Alvarez and Gaioli, 1998). On the other hand, the 
equation of motion corresponding to (Nq) is a generalized- local-in-time form of the master 
equation, with time-dependent coefficients (this was anticipated in I but finally showed in 
Garcia Alvarez and Gaioli, 1998). 

In Section 2 we reobtain the generalized Langevin equation (1.42) but in this case without 
considering it in mean values. Then a new term appears, which plays the role of a stochastic 
acceleration. Evaluating the time-dependent coefficients through a perturbative procedure 
up to the first relevant order we recover the standard form of the Langevin equation, where 
the coefficients are independent of time and the stochastic acceleration has a zero-centered 
distribution leading to white noise in the classical limit. 

In Section 3 we show that the same perturbative analysis mentioned above performed on 
equation (1.33) leads to the solution of an approximated equation for the mean population 
{Nq,) , introduced by van Kampen in his 1992 book [equation (XVII. 2. 30) of van Kampen, 
1992]. Van Kampen's equation, which is the Born approximation of the generalized master 
equation derived in a previous work (Garcia Alvarez and Gaioli, 1998), makes explicit the 
temporal behavior of (Nq) : It decays exponentially until reaching thermal equilibrium with 
the heat bath, i.e. the Bose distribution. This behavior has been shown in the figures of 
paper I, obtained from the exact solution of the model. However, a carefully evaluation of 
the asymptotic value of (Nq) shows us that some deviations from the Bose population arise. 
That is, at high and intermediate temperatures the equilibrium value corresponds to the 
Bose distribution at the renormalized frequency Q + 5Q, where 5Q is a shift proportional to 
the square of the perturbation strength, while at low temperatures a power-law behavior is 
found. This deviation is related to the long-time tail of the decay probability PQn(t) -known 
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as Khalfin effect (Khalfin, 1957)- of the unstable one-particle state \Q) [see equation (1.34)]. 

In Section 4 we outline our main conclusions and there are also two appendices: Appendix 
A contains the details of the calculation of the autocorrelation function of the stochastic 
acceleration and Appendix B includes an explicit derivation of the exponential decay law 
and the Khalfin effect. 

II. GENERALIZED LANGEVIN EQUATION AND THE "STOCHASTIC" 

ACCELERATION 



Let us remember that the Hamiltonian [equation (1.5)] of the composed system is given 



by 



/ l\ n / In n 
H = Q [B^B + - + ]>> n (bib n + - + £ g n (Bb{ + B%) 

V Z " / n=l V Z " / n=l 



The position of the subsystem Brownian oscillator X(t) = (2MQ) 2 B\t) + B(t) 
rewritten, from equation (1.26), as 



(1) 



can be 



X(t) = a{t)X(0) + b(t) 



MVL 



/(*), 



where 



N 

a(t) = \®u\ 2 cos(a u t), 

N 

b(t) = J2 \& v \ 2 sm(a v t), 

1 9r , 
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n=1 ^a u -iu n 
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Since we have two constants of integration X(0) and -P(O), and a particular solution f(t), 
X(t) satisfies a second order differential equation such as 



X (t) + Q 2 (t)X(t) + r(t) X (t) = F{t), (3) 
with the inhomogeneous term given by 

F(t)=f(t) + tf(t)f(t) + r(t)f(t). (4) 

The unknown coefficients fi 2 (t) and T(t) can be easily determined solving the linear system 
which results by replacing the two independent solutions of the homogeneous equation: 

a (t) + n 2 (t)a(t) + r(t) a (t) = 0, 
b(t) + Q 2 (t)b(t) + T(t)b(t) = 0, 

that is 



ab —o a ab —o a 

Garcia Alvarez and Gaioli (1998) have written these coefficients in terms of the survival 
amplitude (and its complex conjugate) of the state \u ) = \Q) = |0): 

Ann(t) = (n \e~ iht \ n) = a(t) + ib(t), 

N N 

where h — (Q\ + w„ \u n ) (uj n \ + J2 9n (|^) {w n \ + \w n ) + C [equation (1.7), we 

n=l n=l 

are calling now the Hamiltonian h instead of Hi, according to the notation of the Garcia 
Alvarez and Gaioli (1998) work], and C is the zero-point energy. This survival amplitude is 
the cornerstone of the theory of unstable quantum systems. 

On the other hand, the mean value of the "stochastic" acceleration F(t) vanishes 
since in the thermal initial distribution (1.29) [p(0) = ps(0) ® — e * -" b H ,\ , where H b = 
Y^n=i u n (b\,b n + 1/2) and tr& is the partial trace over the reservoir], the bath operators 
have vanishing mean values, 



(F(t)) = 0. (6) 
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Performing a perturbative expansion up to the first relevant order in the coupling pa- 
rameter we can recover the standard form of the Langevin equation. So, in this case we 
have 



x{t) + {tt + 5nyx{t) + 1 x{t) = F{t), (7) 

where Sfl is the shift of the frequency and 7 the damping coefficient. 

To see this we evaluate the survival amplitude up to the second order from (Sakurai, 
1995) 

Ann(t) = e~ int cnn, (8) 

where 

n r t r t' 

can = 1 - £ 9k / / e^'e^'dt'dt", (9) 

and ujQ k = Q — u k . Derivating equation (g) with respect to time, up to second order, we 
obtain 

N 



l-^Y,9t (-» I e^ T dr 



Performing a long-time approximation (t > ^) 

lim (-% [ e^dr) = 6+(a) = — = PP- - in5(a), (10) 
t-«o \ Jo J a + ie a K ' v ' 

where PP denotes the principal part. Therefore the survival amplitude of the state has 
a decaying exponential contribution [cf. equation ( |B10|) 1 



A nn (t)=e-i n+m -^y, (11) 

whereQ 



2 The following expressions must be understood as if a continuous limit has been taken in such a 
way that summations become integrals. On the contrary, we must replace the principal part and 
delta distributions by their corresponding approximants. 



N 




5Q = PP 



(12) 



and 



N 

7 = 27^ -fi). (13) 

If we now use the second order value (|TT|) of the amplitude Ann into equation (^), a 
straightforward evaluation leads to 



so we retrieve standard expressions usually derived from the Born and Markovian approx- 
imations. In our case the Markovian approximation is not necessary because equation @ 
is just local in time, in contrast with the more known integro-differential form commonly 
found in the literature (Louisell, 1973; 1977; Sargent et at, 1974; Lindenberg and West, 
1984; 1990; Ford et al, 1988; Meystre and Sargent, 1991; Cohen- Tannoudji et ai, 1992; 
Mandel and Wolf, 1995). 

The autocorrelation function of the "stochastic" acceleration is defined by 



Up to the second order, by means of a long but straightforward calculation (see Appendix 
A), we obtain 



n(t) = Q + SQ, 



r(t) = 7, 



(14) 



K(t) = ±(F(0)F(t) + F(t)F(0)). 




m=l 



(15) 



For an initial thermal distribution for the bath oscillators, i.e. 2 (N m (0)) + 1 



coth^, 



and by taking the limit of a continuous bath (see paper I), K(t) becomes 




where 



g 2 {w) is defined as g 2 (u)Auj 




g 2 (Ullersma, 1966; van Kampen, 



1992). 
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Considering a coupling function g 2 (u) which is peakedQ around uo = Q, then we can 
replace |g 2 (co>) by ^27rg 2 (Q) in the integrand. Taking into account that, up to the 

second order [see equation ([Op], the damping factor in the Langevin equation is 7 = 2ng 2 (Q), 
we finally have 

7 f°° , * .1 Pi™ . lk B T d ^ fnk B Tt\ 



*Mfr) ~ tt- / acuAkucoth cos cut = — — — - coth 

w Mix Jo 2 M dt 



In the classical limit (h — > 0) it goes to the classical stochastic, delta-correlated, autocorre- 
lation function 

L m o^) = ^w- < 16 > 

originally proposed by Langevin. 

Equation (|16|) corresponds to instantaneous correlated fluctuations, which leads to a 
Markovian process (instantaneous memory loss, since the values of the stochastic acceleration 
at two different times are not correlated). F(t) is the source of noise (fluctuations), known 
as white noise. In this limit the usual Langevin equation ([/]) together with properties (^) 
and ( IB ) are recovered from a general formulation. A new ingredient is that in equation (|7J) 



the oscillator frequency is shifted in a value Sfl. 

III. BEHAVIOR OF THE MEAN POPULATION OF THE BROWNIAN 

OSCILLATOR 

Van Kampen (1992) has shown that the elimination of fast microscopic variables leads to 
a closed expression for the mean value of the occupation number of the Brownian oscillator 

f t (Nn(t)) = -^(Na(t)) + ^- [ , (17) 



3 This condition is necessary in order the rotating wave approximation to be valid (see, e.g., Gaioli, 
1997). 
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where irrelevant variables get in only through the initial thermal state at a temperature 
T = l/k B p. 

In what follows we use the time-dependent perturbation theory to derive equation ( |TFD 
starting from equation (1.33) of I: 

N 

(N n (t)) = Pnn(t) (N n {0)) + £ Pn n (t) (N n (0)) . (18) 

n=l 

We consider that the perturbation [interaction of equation ([I])] is time independent and 
turns on at t = 0. In this case the Dyson series for the transition amplitude to second order 
can be written as (Sakurai, 1995) 

A nm (t) = e m {c[^ n + c[^ n + c^ n + . . .J , (19) 

where 

c (0) = S 



<£L = -i fe^'vmndt', (20) 







•v r t rf 



k=0 Jo Jo 



with v mn = (u> m \v\v n ) = SVrA.o + 9nSm,o, except for vqo = 0, and u mn = oj m - uj n . The 
evaluation of the transition probabilities to second order gives 

with 

r nm = 2irvl m 6 t (u n -(jj m ), forn^m, (22) 

and 

T nn = -2tt v lm<>t(u n - uj m ), (23) 
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where 8 t (a) = %f is a function which approaches Dirac's delta when time goes to infinity.^ 
Introducing these results into equation ( |l8l) and taking into account that r o = r^n = —7 
and that, in this case, Pqq ~ 1 — jt, we obtain 

(Nn(t)) ~ (1 - jt) (Nn(0)) + jt (iV«^=n(0)) 

= (1-^X^(0)) + ^^, 

which is the second order perturbative expansion of 

(tf n (f)> * e"^ (iV n (0)> + (l - e-*) — ^ (24) 

Equation (|24j) corresponds to the solution of equation (|T7| ) [cf. equation (XVII. 2. 29) of van 
Kampen (1992)]. This is also the second order approximation of the solution of the exact 
master equation of Garcia Alvarez and Gaioli (1998). 

We have to emphasize that this is a perturbative result. A full exact calculation can 
be carried out in the limit of a continuous bath (as it was taken in Section 6 of paper I). 
We now calculate the asymptotic long time equilibrium value of the mean population of the 
Brownian particle under this limit. We begin with equation (1.76) but setting u m i n and uj ma , x 
equal to zero and infinity respectively, so we have 

(JV„(oo)> = T du ^ 2 1 ■ (25) 

By considering now "no-low" temperatures and a small coupling it is easy to see that 
g 2 {ui) Rl 1 ^) picks the value at u = Q + 5fl and therefore has a 5{uj — Q — 5Q) behavior. 



4 Let 77 be the width of the interaction function g^. Since the function 6t (a) has a width 4ir/t, in 
order to behave as a delta distribution we need times such that t ^> An/rj. This is the meaning of 
long times. 

We also see that the times involved in this approximation must satisfy t <C I/7, where 7 = —Too, 
in order equation (^Tj) to be valid. This condition and the long-time limit restrict the time-range 
to 47r/7/ <t < I/7. 
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Thus equation (p5]) leads to the thermal equilibrium value at the shifted frequency Q + Sfl, 
i.e. 

(iVn(oo)) = (N„ =Q+SQ (0)) = e ^ (n+ ^ ) _ 1 - (26) 

In the classical limit (iVn(oo)) ~ [(3(Q + SQ)]^ 1 , which leads to the equipartition of 
energy E « ksT (two quadratic degrees of freedom in a one- dimensional configuration 
space) and the heat capacity C v ~ ks- 

At low temperatures, deviations (an inverse power-law falloff in expectation values) from 
this equilibrium distribution were already reported by some authors (Lindenberg and West, 
1984; 1990; Haake and Reibold, 1985; Joichi et al, 1997), since the effect of the coupling 
becomes macroscopically observable. This effect is similar to the deviation from the expo- 
nential decay law for long times, which was first described by Khalfin (1957). In Appendix 
B we explicitly calculate the survival amplitude using the resolvent method. In such a case 
we show that Aqq reduces to [equation (|B8"D 1 



Aoa{t) = r duJ g2{U) 2 e ~ lWt - ( 2? ) 

|#;V)| 

Comparing equation ( p7|) with the low-temperature limit of equation (|25|) , 



(iV n (oo)>= [°°du ^ 2 e~^, (28) 







we can see that, by making the identification (3 = it, both expressions are equivalent. 
That is, the low-temperature regime of the mean occupation number corresponds to the 
long-time behavior of the survival amplitude of the unstable state \Q) . This is one of the 
deep interrelationships between statistical properties of nonequilibrium ensembles and the 
behavior of unstable quantum systems. This was possible since our model [equation ([[])] can 
be decomposed by sectors of fixed number of quanta [see equation (1.6)], which is one of the 
reasons for our choice of the model. 

We then use the long-time limit of the survival amplitude (Khalfin effect) of Appendix B 
[equation ( |B12[ )] in order to estimate the low-temperature anomalous behavior of the mean 
population, namely 
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(JVh(oo)) ~ (k B T) n+1 . (29) 

The energy of the Brownian oscillator behaves like E ~ T™ +1 and so C v ~ T n , which 
is an analogous result to the low-temperature behavior of the Debye model of phonons 
(Huang, 1987). One of the advantages of this kind of calculation is that equation (|29| ) can 
be experimentally measured and then it provides an indirect proof of the existence of the 
Khalfin effect, which is very difficult to measure because of the time scale involved. However, 
the other deviation from the exponential decay law, the Zeno effect (Misra and Sudarshan, 
1977), was recently measured for the first time (Wilkinson et al., 1997). 

We have seen that the anomalous behavior of the mean population at low temperatures 
[equation (|29|) 1 is related with the deviations of the exponential decay law of the unstable 
initially prepared state \Q) at vey long times. We can also see another suggesting relation 
between this anomaly and a generalized statistics proposed by Tsallis (1988) ten years ago 
(see also Curado and Tsallis, 1991). Let us see the origin of this conjecture. 

Biiyiikkilig and Demirhan (1993) and Biiyiikkilig et al. (1995) have shown that for a 
set of bosons, labeled by the index k, the mean occupation number corresponding to the 
generalized Bose-Einstein canonical distribution is approximately^ given by 

where q is a parameter which characterizes the nonextensive nature of the system, and thus 
depends on the long-range nature of interactions present in the system. This parameter 
is such that for q — > 1 one retrieves the standard results, i.e. the Bose-Einstein mean 
population. For q ^ 1, in the low-temperature regime, we can neglect the terms independent 
of temperature in the denominator of expression (|30|). Then, we have 

(N k ) ~ {k B Tf^ . (31) 



'The fact that this is a non-exact result was noticed by Pennini et al. (1995). 
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Compare equation ([31] ) with equation fl2~9|) . We see that the equilibrium distribution 
reached by the Brownian oscillator resembles that of the bosons in thermal equilibrium for 
a nonextensive system, according to Tsallis' prescription. From these equations we obtain 

n + 2 

q = — TT' 
n + 1 

a number which satisfies 1 < q < 2, since n > 0. Maybe the explanation of this behavior 
is that the Brownian oscillator is not able to cover all accessible quantum states (as a 
consequence of strong quantum correlations at low temperatures) and then it cannot reach 
the most probable distribution according to the ergodic hypothesis. 

IV. CONCLUSIONS 

Along this paper the autocorrelation function of the stochastic acceleration and the 
asymptotic mean population of the Brownian oscillator were analytically evaluated from a 
deterministic quantum dynamics. As regards the Langevin equation we have provided the 
stochastic term which was skipped in paper I. For the mean occupation number we have 
found that it reaches thermal equilibrium at the bath temperature, corresponding to the 
Bose population. At low temperatures a deviation from this population was found, which 
has a common origin with the deviations from the exponential decay law. However, the 
Khalfin effect is very difficult to measure since usual observation times of unstable quantum 
systems are much shorter than the time the decay law is no longer exponential. 
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APPENDIX A: AUTOCORRELATION FUNCTION OF THE STOCHASTIC 

ACCELERATION 



According to equation K(t) is given by 

K(t) = \ [(/ (o) / (*)+ / m 2 (t)f(t)+ / (0)r(t) / (f) (Al) 
+fi 2 (o)/(o) / (t) + n 2 (o)f(o)n 2 (t)f(t) + fi 2 (o)/(o)r(t) / (t) 
+r(o) / (o) / (t) + r(o) / m 2 (t)f(t) + r(o) / (o)r(t) / (t) 
+ (o ^ t}} , 

where (0 t) stands for interchanging t = with £. /(£) can be written in terms of the 

1 N 

/(*) = £ KWC(O) + /i-c] • (A2) 

Considering second order contributions only and taking into account that Sfl and T are 
time independent up to this order, we can rewrite equation ( |A"I| ) as 

K{t) = \ [(/ (o) / (t) + n 2 f (o)f(t) + n 2 f(o) / (t) + fi 4 /(o)/(t)) 
+ (o ^ t>] , 

since /'s are linear in A^m's [see equation (|A^)]. We need to expand the amplitudes up to 
the first order. From equations (|19D and (^) we have 

Jo 

Solving the integral it is 

uj m - i I \ ' 
The second derivative of Afi m (t), appearing in /, is given by 

V m n ( 2 -iuj m t 
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Taking into account that An m (0) = 0, from which /(0) = 0, then we must only calculate 



K(t) = I \(f (o) / (t) + n"f (o)/(t)) + (o <- t) 



Let us see each term step by step. Let K\ and K 2 be defined by 



Ki{t) = -(f (0)f(t)+f(t)f(0) 



and 



^ 2 (t) = y Q (°)/(*) + /(*) / (°) 



Therefore 



11 ^ 

Kl{t) = 22Mti ^ (°) 4 ^ ( & L^') (°)+ 4 (°) ^' (*) (^4) (°) 



m,m'=l 



+ An m (t) A* nm , (0) (blb m ) (0)+ (*) An m < (0) (& m &U (0) 



which, from equation (1.30) and the commutation relations (1.4), is reduced to 
1 N 

Ki(t) = ^ £ {Re [An™ (0) (t)\ [2 (iV TO ) (0) + 1 
1 



m=l 
N 



N ( II ~l 

£ [2 (N m ) (0) + 1] (u m + n) U m C osu m t - n 2 cosftt) 



An analogous calculation for K2(t) leads to 

Q AT 

= ^7 E { Re pin™ (0)^ m (t)] [2 (iV m ) (0) + 1]} 



2M 
1 

2MVL 



m=l 

N 



]T <j [2 (JV m ) (0) + 1] \ Vmnl ( Um + fi) (q 2 cos fit - ft 2 coscj m t 



m=l 



Joining K\ and K 2 we finally obtain equation 



1 N 

K(t) = {[2 (N m ) (0) + 1] \v mn \ 2 (uj m + ft) 2 coscu m t} . 



m=l 
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APPENDIX B: EVALUATION OF THE SURVIVAL AND TRANSITION 

AMPLITUDES 



We analyze the analytical structure of Ann(t) in order to show how an exponential 
contribution arises for a significant range of time and how deviations from this behavior 
appear. Using the well-known identity between distributions 

' PV-^i7r5(x), (Bl) 



x ± ie x 

we can obtain, for x = u — h, the following integral representation of the evolution operator 

1 r°° 

-iht _ I j, ,-iwt 



e - tM = / du)e -^t [G_(u) - G+(u)] , (B2) 



where G±(a) = a ±l e _ h is the retarded (+) [advanced (— )] Green function of the time inde- 
pendent Schrodinger equation (the resolvent), corresponding to the total Hamiltonian. From 
( |B2| ) we can evaluate the survival and transition amplitudes. Such calculation involves the 
knowledge of the partial resolvents (Schwinger, 1961; Messiah, 1962), departing from 

(a±ie-h)G±{a)=I. (B3) 

By taking the matrix elements in equation ( |B3| ) we have 

roo 

(a±ie- fi) (fil G±(a) Ifi) + / dujgiu) U G±{a) \Q) = 1, (B4) 

Jo 

(a ±ie - u) {u\ G±(a) \Q) = g(u) G±(a) |fi> = 0. (B5) 
So, from ( |B4|) and (|B5|), the desired matrix elements of the resolvent are 
(tt\ G±{a) \Q) = ztt- = R±(a) 

J U a±ie 



and 

(w| G±(a) |0) = f u) (0| G ± (a) |fl> . 
a ± ie — uj 

Now returning to (|B2|) we can obtain the survival and transition amplitudes as 
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A nn (t) = (Q\e~ iht \n) = ^- 



doj'e-^ [R-(oj r ) - R+(oj') 



(lu\ e 



-iht 



duj e 



I — iia't 



g{u')R_{u') g(u')R+(u') 



to — it — UJ 



u> + it — uj' 



(B6) 



(B7) 



i.e. taking the Fourier transform of the difference between the advanced and retarded reduced 
resolvents (note that we have not introduced any initial condition for the amplitudes). It 
can be proved that the reduced resolvents have the general expression 

R ± (uj) ' 



uj ± it — Vt — X±(u;) 

where S(a;) is the level-shift operator in the subspace generated by \Q) . The reduced resol- 
vent R+(uj) [R-(uj)} is the analogous of the exact Feynman (Dyson) electron propagator with 
Ti±(uj) playing the role of the Dyson (1949) mass operator. In our case £±(a) = J °° duj ( 9 
Using ( |B1| ) we can rewrite it as £±(a) = A(ct) =f with 



aiie—LU ' 



Ma) 



2 ' 

PV fcfc/M, 
Jo a — uj 



7(a) = 2-ng (a), 

which are nothing else than equations flT2"| ) and ( |T3D in the case of a continuous bath. Taking 
into account that g{yj) = for uj < 0, equation ( |B6| ) can be rewritten as 



duj- 



-iu)t 



B8) 



R+ x {uj 

In the theory of unstable states (Messiah, 1962; Goldberger and Watson, 1964; Cohen 
Tannoudji et al, 1992) it is common to find Aqq written as 



A m (t) 



1 

27 



duj- 



7(w) 



-iuit 



(B9) 



[uj - n - A(u;)]' + | 7 » 

which allows one to easily study the different decay regimes. If j(oj) is small, the term 
[uj — Vt — A(uj)] 2 is large compared with 7 2 (cu) except when uj ~ Q + A(O). Thus we replace 
7(0;) by 7 (Q) = 7 and A(uj) by A(f2) = 5Q. The Lorentzian function resulting from this 
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replacement is known as Breit-Wigner (1936) distribution. In this case equation (p9| ) has 
an analytical result 

Ann(t) = e -W^ e -i', (BIO) 

which, as it is expected, retrieves the well-known exponential decay law, as it was originally 
derived by Weisskopf and Wigner (1930). Some deviations from this exponential decay arise 
as we inspect equation ( |B9| ) more carefully. If we retain j(oo) in the numerator of (|B~9|), we 
can rewrite this equation as 

A nn (t) = — [ du ZM e -^. (Bll) 

W 27r J / j [ u - n - 5tt] 2 + ± 7 2 7 

The r.h.s. of equation (|B11| ) is the convolution product of ( |B10|) and the Fourier transform 
of 7 _1 7(c<j). Since 7 (cu) has a finite width, its Fourier transform 7 (t) has also a finite width. 
We also have that 7 (co>) is null for u < and is not infinitely differentiable at uj = 0. Then, 
if we suppose that 7 (co>) goes as uj n for small u, then j(t) behaves like t~( n+1 ) for very long 
times. This is known as Khalfin (1957) effect, namely 

Ann(t) ~ (B12) 
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